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Abstract 

The software alignment of planar tracking detectors using samples of charged 
particle trajectories may lead to global detector distortions that affect ver- 
tex and momentum resolution. We present an alignment procedure that 
constrains such distortions by making use of samples of decay vertices re- 
constructed from two or more trajectories and putting constraints on their 
invariant mass. We illustrate the method by using a sample of invariant-mass 
constrained vertices from D° — » K~n + decays to remove a curvature bias in 
the LHCb spectrometer. 
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1. Introduction 

The calibration of the position and orientation of tracking detectors in 
high energy physics experiments is called alignment. The input for alignment 
comes from two sources, namely survey information collected during assem- 
bly or after installation, and hit residuals of reconstructed charged particle 
trajectories (tracks). With track-based algorithms an alignment accuracy 
can be reached that well exceeds the single-hit resolution. 

The track-based alignment algorithms considered here optimize the total 
fit quality of a sample of tracks, for example the total track fit x 2 , with respect 
to a set of numbers that parametrize the misalignments. The parameters are 
usually chosen to be the positions and orientations of individual detector 
elements. We denote the set of alignment parameters with a generic symbol 
a. The condition that the total x 2 be minimal with respect to a can then be 
written as 

tracks i 

where the sum runs over all tracks in the calibration sample. 

A common problem in the application of track-based alignment algo- 
rithms is related to so-called weak modes. These denote alignment degrees 
of freedom to which the total track x 2 is mostly or completely insensitive. A 
global translation or rotation of all detector elements with respect to a com- 
mon point (i.e. a coordinate transformation) is a perfect weak mode. Less 
trivial weak modes are related to global distortions and depend on the de- 
tector geometry. In parallel plane detectors, shearings, such as that depicted 
in Fig. [TJ are also examples of weak modes. 




Misaligned position 

Figure 1: Illustration of a shearing in a planar detector. 

One reason weak modes are a concern in alignment is that they lead to 
poor convergence. An approach to treat weak modes is to impose additional 
constraints. For example, information on alignment parameters obtained 
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from detector survey can be used in the minimization procedure by including 
additional penalty terms in the x 2 . 

Weak modes related to global distortions pose a particular concern be- 
cause they can lead to biases in track parameters that affect the performance 
of an experiment. A global translation of the entire detection apparatus 
changes the numerical values of the track parameters without changing the 
kinematics relevant for a physics analysis, such as an invariant mass or a 
decay angle. However, global distortions that affect the relative position or 
direction of tracks will introduce a bias in kinematic observables and degrade 
the overall detector resolution. 

A particularly interesting weak mode related to a global distortion is 
the so-called curvature bias. This weak mode appears both in cylindrical 
detectors with a solenoidal magnetic field (in which it is sometimes called 
sagitta bias or curl) and in forward spectrometers with a dipole magnet, 
as illustrated in Fig. [2] The momentum of a charged particle is measured 
via its curvature radius, which requires at least three measured coordinates. 
In a uniform magnetic field B, the inverse of the curvature radius, which 
we shall call the curvature u, is related to the momentum component p± 
perpendicular to the field by 



where Q is the charge of the particle. The weak modes illustrated in Fig. [2] 
introduce a bias in the curvature 



that is approximately independent of particle momentum and direction. A 
constant curvature bias leads to a momentum bias that depends both on 
momentum and charge. As we shall see later, one signature of a curvature 
bias is a shift in the reconstructed invariant mass of a two-body decay that 
is proportional to the momentum difference between the two final-state par- 
ticles. 

The weak mode that leads to a curvature bias only exists in the presence of 
a the magnetic field. Consequently, a curvature bias can be constrained with 
field-off data, provided that detectors do not move if the field is turned on. 
Unfortunately, the latter condition is not often fulfilled. In central detectors 
a curvature bias can be constrained with samples of reconstructed cosmic ray 
tracks by exploiting that the 'top' and 'bottom' leg of the track must have 
the same curvature [I]. In forward detectors, such as the LHCb detector at 
CERN [2], this technique does not work and alternative methods must be 
deployed. 

We report here on a novel method to constrain the curvature bias by in- 
cluding a x 2 contribution from a mass-constrained vertex fit of a multi-body 




(3) 
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Figure 2: Illustration of a curvature bias in a cylindrical detector geometry (left) and in a 
forward detector geometry (right). 

decay. Implementations of kinematic and vertex constraints for alignment 
have been presented before [U El Hj. These applications rely on parametriza- 
tions of a two-body decay that cannot easily be extended to more than two 
tracks. Our method does not need a special track parametrization or track 
fit and can be used with vertices with any number of tracks. 

The outline of this paper is as follows. In Section [2] we briefly discuss 
the track-based alignment procedure, referring the reader to a previous pub- 
lication for details. In Section [3] we present the implementation of a vertex 
fit for use in the alignment. In Section [4] we discuss the use of vertex con- 
straints in the alignment of the LHCb detector. As an illustration we compare 
the performance of alignments obtained with and without mass-contrained 
D° -> K-TT+ decaysQ 

2. Alignment with tracks and vertices 

The application of a x 2 minimization procedure for track-based align- 
ment is extensively discussed in the literature (see e.g. Ref. [5] and references 
therein). For our discussion we use the formalism presented in a previous 
publication [6], in which we have shown how to use a Kalman-filter track 
fit in a closed-form alignment procedure. We follow the notation from that 
paper and from a seminal paper on the application of Kalman filters in high 
energy physics [7j. 

We denote the set of alignment parameters by a and the track parameters 
of track i by X{. For a track with N hits, we write the N dimensional vector 
of hit residuals schematically as 




(4) 



throughout this paper charge-conjugated modes are implied. 
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where is the vector of hit coordinates, each of which is a function of 
(a subset of) the alignment parameters. The function hi is often called the 
measurement model. It expresses the expected hit coordinates in terms of the 
track parameters^] The x 2 of track i can now be written in matrix notation 

as 

X- = rf V' 1 Tl , (5) 

where V{ is the NxN covariance matrix of the measurement coordinates. The 
latter is usually diagonal. The best-fitting track parameters are obtained by 
minimizing xf with respect to Xj for a given value a of the alignment param- 
eters. Since the measurements depend on a, so do the track parameters. 

While the track parameters are different for each track, the alignment 
parameters are common. To obtain the optimal alignment one minimizes 
the total x 2 °f a sample of tracks simultaneously with respect to the track 
parameters X{ and the alignment parameters a. By taking into account how 
the track parameters depend on the alignment parameters, the problem can 
be reduced to a minimization problem with the dimension of a only jSJ El QUI 

EJ. 

Starting from an initial alignment ao, the solution for a = ao + Aa is 
obtained by solving the set of linear equations 

dV Aa = - ^ . (6) 



da 2 



da 

a 



The first and second derivatives are obtained by summing the contributions 
from all tracks and can be expressed as [101 16] 

tracks i 

and 

tracks i 

where Ri is the covariance matrix of the residuals after the track fit. The 
latter is expressed as 

Ri = Vi - HidHT , (9) 

where Hi is the derivative of with respect to the track parameters of track 
i and Cf is the covariance matrix for the track parameters. To obtain the 
expression in Eq. \f\ one exploits the fact that the x 2 contribution for each 



Note that one can choose where to put the alignment parameters a in Eq.|4| The 
residual can also be written as rj = to.; — hi(xi,a) [6]. Only the residual and the 
derivatives of the residual to track and alignment parameters enter the formalism. 
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track has been minimized with respect to the track parameters for the initial 
alignment a [6]. 

One ingredient to Eqs. [7] and [8] are the derivatives of the residuals (or 
measured hit coordinates) to alignment parameters]^] Their computation 
depends on the implementation of the detector geometry and is outside the 
scope of this paperQ Another ingredient is the track parameter covariance 
matrix d. In LHCb we use a Kalman filter track fit that takes multiple 
scattering and energy loss into account and follow the approach derived in 
Ref. [B] to compute C$. 

In Section 4 of Ref. j6] the formalism above is extended with vertex con- 
straints. The proposed method relies on a vertex fit. The vertex fit computes 
new track parameters on the assumption that all tracks in the fit originate 
at a common point. Using the covariance matrix of the tracks the difference 
between the track parameters before and after the vertex fit can be propa- 
gated to the hit residuals r,. The vertex fit introduces a correlation between 
the parameters of different tracks. Therefore, the x 2 contributions from dif- 
ferent tracks are no longer independent. This means that the residual vector 
rj now spans residuals from all tracks included in the vertex and the covari- 
ance matrix Cj is now the covariance matrix for the parameters of all tracks 
included in the vertex. The method is valid for any number of tracks in a 
vertex. The strength and novel aspect of the method is that no special track 
fit is required. The one ingredient that is missing in the discussion in Ref. [6] 
is the vertex fit itself. We present the vertex fit in the next section. 

The actual alignment procedure now consists of the following steps. First, 
the tracks are fitted. Next, subsets of tracks that are identified to come 
from a single vertex are combined with a vertex fit. Using the formalism 
in Ref. [6] a single residual vector r; and corresponding covariance matrix 
Ri are computed for each multi-track object and added to the derivatives in 
Eqs. [7] and [8] Selected tracks that are not used in a vertex can also be added 
to these derivatives. Finally, once all contributions in the sample have been 
accumulated, new alignment parameters are computed using Eq. [6] 

3. The vertex fit 

A vertex fit combines the trajectories of a set of charged particles with the 
constraint that the particles originate from a common point. The input to the 
vertex fit are the reconstructed parameters and covariance matrix for each 



3 We denote these as partial derivatives since at this stage one ignores the contribution 
to the derivative that comes through the track parameters. 

4 Their evaluation in LHCb is similar to that discussed in Section 3 of Ref. [TU] for the 
ATLAS detector. 
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of the tracks. The output of the vertex fit is a vertex position, a momentum 
vector for each of the tracks and a corresponding covariance matrix. 

For the implementation of the vertex fit we use the Billoir-Friihwirth- 
Regler algorithm [TTJ. We follow the notation used in Ref. [7J, with small 
modifications to remain consistent with the symbols used for the track pa- 
rameters in the previous section. For the application in alignment we extend 
the formalism with a mass constraint. 

Tracks are locally parametrized by a 5-D vector, generically denoted by 
the symbol Xi where the label i enumerates the tracks in the vertex. The 
covariance matrix of the track is denoted by Cf and its inverse by Gi. We 
denote the 3-D vertex position vector with the symbol £ and the 3-D vector 
that parametrizes the momentum vector of outgoing track % with q^ The 
measurement model, qi) expresses the parameters of track i in terms of 
£ and qi. The residual of track % is then defined as 

n = Xi - hi(£, q { ) . (10) 

Note that the symbol Xi used for the track parameters is the same symbol 
that appears in Eq. [4] and beyond. However, the role of the track param- 
eters is different: in the track fit, the track parameters are free parameters 
determined by the track fit. In the vertex fit, the track parameters and their 
covariance matrix are input to the fit. In the following we show how the 
track parameters are changed if a vertex constraint is applied. 

The x 2 °f th e vertex fit is written as 

X a = X>?G,r,, (11) 

i 

where the sum runs over all tracks in the vertex. The solution to the vertex 
fit is the set of parameters (£, qi, . . . , q^) that minimizes this x 2 ■ 

As in Ref. [7J [TT] we linearize the measurement model around the current 
estimate (£o,<?i,o), 

hi{tqi) = hi{i Q) q ifl ) + AiH-io) + Bi{ qi -q ifi ). (12) 

The measurement model and its derivatives Ai and Bi follow from the para- 



metrization chosen for tracks and vertices. In Appendix A we shall discuss 
a definition suitable for a forward detector as LHCb. The total derivatives 
of the x 2 with respect to the vertex position can now be written as [11] 

^ = -2XXGfr, and % = ^^GfA ll (13) 
where we introduced 

Gf = Gi-GiBiWiBjGi (14) 
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with 



The updated vertex parameters are obtained with 



dV\ 1 dx 2 



° V d£ 2 J d£ 
while the covariance matrix for £ is given by 

'dVV 1 



(15) 



(16) 



Cov(0 = C c = 2 



(17) 



The momentum parameters of the outgoing tracks can be computed with 

Qi = Qifl + WiBfGi(xi - hi(^ q ifi )) . (18) 



The covariance matrix for two outgoing tracks i and j is given by 
Cov(g i; qj) = D id = 5^ + (W^fG^) C^AjGjBjWj) , 



(19) 



where 5^ is one if i — j and zero otherwise. The covariance matrix for the 
vertex position and the momentum of track i is given by 



CovO^)^ = -WiBfGiAiCz. 



(20) 



The vertex fit can be iterated until a certain convergence criterion is met, 
for example a sufficiently small change in the vertex x 2 °f Eq. 11 Note 



that the computation of the covariance matrices in Eqs. 19 and 20 is CPU 
intensive, but not needed for obtaining the x 2 or the position and momentum 
updates in Eqs. 16 and 18 [TTJ. Consequently, their computation can be 



delayed until the fit has converged. 

An invariant mass constraint adds a new term to the total vertex x 2 i n 
Eq. 11 In contrast to the x 2 contributions from the individual tracks, the 
mass x 2 term has non-zero derivatives to all momentum vectors gj, instead 
of just one of them. As a result the optimized formulation of the vertex 
fit above cannot easily be extended with a mass constraint. Therefore, we 
have chosen to add the mass constraint after the fit to the vertex position, 
ignoring small effects due to the non-linearity of the fit. 

For each final state particle i, we compute a relativistic four-momentum 
vector ki, which is a function of the momentum parameters g« and a mass 
hypothesis mj. We denote the (4 x 3) matrix for the derivative of ki with 
respect to by 

K = (21) 
dqi 
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The total four-momentum of all tracks assigned to the vertex and its covari- 
ance matrix are given by 

fctot = }] fei(gi) mi) and Cov(& to t) = / ] Kj A,j Kj, (22) 



where and are the result of the unconstrained fit (Eqs. 18 and 19 



respectively). From k to t we compute the mass and form a residual for the 
mass constraint, 

r M = m(fctot) - m , (23) 

where tjiq is the known mass of the reconstructed decay. Defining the 1x4 
derivative matrix 

the variance of the constraint is given by 

R M = H M Cov(k tot )H^. (25) 

In case the natural width T of the decaying particle is not small compared 
to the invariant mass resolution, such as for the decay Z° — > one can 

add Tq to the variance of the constraint. 

Using well-known expressions for the Kalman filter [7] , the mass-constrained 
vertex position and track parameters now become 



f= E T k K T k HlR^r M , (26) 

k 

q t = qi-J2 Di, k KlH^R^r u ■ (27) 



k 

The updated covariance matrices are given by 

Q = Q ~ ^2 E k K I HmRm HmKiEi , {21 



k.i 



Dij = D i:j - J2 D^K^H^R^HuKiDu , (29) 

k,l 

Ei = Ei — 2j D i k K k H^R^l H M KiEi . (30) 



k,l 



Finally, for application in the alignment, the vertex-constrained track 
parameters and their covariance matrices must be computed. These follow 
from the measurement model as 

%i = h(£,qi) (31) 
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Figure 3: Sketch of the LHCb spectrometer. The proton-proton collision point is located 
inside the vertex detector on the left (z ~ 0). The main component of the field of the dipole 
magnet is parallel to the y axis. The curved lines represent trajectories of a positively and 
negatively charged particle. 



and 

Cm{x u x j ) = AiC c Aj + AEjBf + BiEiAj + BiD^Bj . (32) 

These equations provide the input to Eqs. 18, 19 and 21 in Ref. [6|j^] In case a 
mass constraint is used we replace the vertex position £, the track momentum 
parameters qi and their corresponding covariance matrices with their mass- 



constrained counterparts in Eqs. 26 -30 This concludes the algebra of the 



vertex fit for use in a track-based alignment algorithm. 

4. Application to the alignment of the LHCb spectrometer 

The tracking system of the LHCb detector is an example of a planar de- 
tector with a forward geometry. It is schematically depicted in Fig. [3] and 
discussed in detail in Ref. [2] . Charged particles produced at the interaction 
point bend in the magnetic field of a dipole magnet with a field integral of 
about 4 Tm. Precision vertexing is provided by a 21-layer silicon strip detec- 
tor located in the field-free region around the interaction point. Four more 
layers of silicon strip detectors just in front of the magnet and another 12 
layers of silicon strip detectors and straw tube chambers behind the magnet 
allow for a precise momentum measurement. 

The coordinate frame in LHCb is defined such that the y axis is pointing 
upwards (parallel to the main component of the dipole field) and the z axis is 
parallel to the beam line with positive z in the direction of the spectrometer. 



5 The symbols Xq\ Cq^ and Cq' j> in Ref. [5] translate in our notation as Cov(x,, 
and Cov(xi,Xj), respectively. 
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The x axis is chosen such that (x, y, z) is a right-handed system. The origin 
is located approximately in the center of the vertex detector and roughly 
corresponds to the average interaction point. 

The implementation of a global minimum x 2 algorithm for the alignment 
of tracking detectors in LHCb has been previously discussed in Ref . [121 EG] • 
Since then we have extended the algorithm to exploit vertex constraints from 
primary vertices and from resonances using the techniques outlined above. 

4-1. Primary vertex constraints 

Primary vertices are important for the alignment of the LHCb vertex 
detector to guarantee an optimal impact parameter and decay time resolu- 
tion [13]. The silicon modules of the vertex detector are assembled in two 
detector 'halves' that are positioned on the positive and negative x side of 
the LHC beam line. The fraction of tracks leaving hits in both halves is 
small. Furthermore, tracks that cross detector planes both in front and be- 
hind the average interaction point (at z — 0) are rare as well. Reconstructed 
primary vertices allow to link detector planes at positive and negative x and 
at positive and negative z. 

We have successfully exploited primary vertices in LHCb alignment using 
the algorithm described above. Although the algorithm in principle allows to 
use vertices with an arbitrary number of tracks, a practical problem occurs 
for high track multiplicity. In LHCb primary vertices often contain tens of 
reconstructed tracks. Reconstructed tracks can have up to 40 hits. For large- 
multiplicity primary vertices the computation of the correlations between all 
hits on all tracks is computationally very demanding. In Ref. [6] it was 
suggested to compute only the correlations between the hits nearest to the 
vertex. However, we have found that this can lead to non-positive definite 
contributions to the second derivative of the x 2 an d therefore is not a viable 
solution. 

A working solution has been obtained by limiting the track multiplicity 
in vertices. For application in alignment we divide a reconstructed primary 
vertex into separate vertices with at most eight tracks each. Tracks are sorted 
such that tracks pointing forward, backward, left and right are distributed 
equally over the different vertices. In our framework the computation time of 
the correlations in these vertices is small compared to the overall reconstruc- 
tion and track fitting time. The loss in statistical power due to the splitting 
of the vertex can be compensated by using more events. 

4-2. Invariant mass constraints 

To constrain weak degrees of freedom in the spectrometer, such as the 
curvature bias discussed in the introduction, mass-constrained vertices from 
D° — > K~tt + , Jfij) — > and Z° — » ji + jj," decays are used. The advantage 



11 



of D° — » K~ir + over other resonances is their large abundance and clean sec- 
ondary vertex signature in LHCb. This allows to select samples of thousands 
of events per hour with practically no background. 

In order to illustrate the effect on the alignment and on momentum mea- 
surements, we show a comparison of two alignment strategies, one obtained 
using approximately 300k selected high-momentum tracks, and another us- 
ing in addition the constraint from 80k mass-constraint D° — > K~ir + vertices 
in the same sample. For this exercise all alignment parameters for the LHCb 
vertex detector were fixed, while all detector elements behind the vertex de- 
tector were allowed to move in the x direction and rotate in the xy plane 
around their center of gravity. To constrain weak modes and ensure conver- 
gence, survey information was used by adding for each alignment parameter 
a penalty term to the total \ 2 constraining it to its survey value within the 
uncertainty assigned to the survey. The alignment process started in both 
cases from alignment parameters obtained with early data p2J IT5] . 



With mass 

constraint 

a Without mass 
constraint 



j i i i i i i i i I 

23456789 10 

Iteration number 

Figure 4: Average \ 2 P er degree of freedom at each iteration of the alignment procedure 
for the alignment with (black line) and without (red squares) D mass constraint. 

In both scenarios the minimization converged quickly as illustrated by 
the average \ 2 per degree of freedom versus iteration, shown in Fig. |4j Al- 
though the average track x 2 after the alignment is practically identical in 
the two cases, the performance in terms of invariant mass resolution is very 
different. Figure [H] shows the invariant mass distribution of D° — > K~ir + and 
J/i/j — > fi + n~ candidates on independent data sets using the two alignment 
sets as input. The mass resolution obtained with the alignment that exploits 
the D° — y K~tt + mass constraint is considerably better. 

To understand this behaviour we consider the effect of a curvature bias 
on the invariant mass of a decay to two oppositely charged particles with 
momenta p- and p + and masses m_ and m + , respectively. Assuming small 
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Figure 5: Invariant mass distribution of D° — > K~ir + (left) and J ftp — > (right) can- 

didates obtained with the alignment using the D° mass constraint (black solid points) and 
the alignment based only on tracks (red open points). 



masses with respect to the momenta, i.e. ignoring terms of order mf/pf, the 
invariant mass is given by 



m = ym 2 + m\ + 2p_p + (l — cos9)/c 2 , (33) 

where 9 is the opening angle between the two particles. As a result of a 
curvature bias (5^, the momentum changes as 

— RH4-w)< (34) 

where we have ignored higher order terms in S w /oj. Note that the sign of the 
bias in the momentum is opposite for the positive and the negative track. 
In a forward spectrometer the momentum component perpendicular to the 
field dominates the total momentum. The change in the invariant mass then 
becomes, to first order in 5 u /oj, 

1 - {p+ - p - ] — 2^ — ± ik) ' (35) 

where e is the positron charge. In other words, one expects a bias in the 
mass that is approximately proportional to the difference of the momenta of 
the two final state particles. 

This effect is demonstrated in Fig. |6j The figure shows the position of 
the peak of the invariant mass distribution for D° — > K~n + and J/ip — > fi + fi~ 
decays as a function of the momentum difference of the final state particles in 
the two scenarios. Indeed, if the mass constraint is not used a clear evidence 
of a curvature bias is observed. The removal of the curvature bias explains 
the difference in mass resolutions shown in Fig. [5] 
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Figure 6: Position of the peak of the invariant mass distribution for D° -> K^tt + (left) 
and J/ifi — > (right) candidates as a function of the momentum difference of the 

two daughter tracks evaluated with the alignment using the D° mass constraint and the 
alignment based only on tracks. 



From the slopes of the graphs in Fig. [5] one can estimate the value of 
b^j^eB. We verify that the two slopes are consistent with a single curva- 
ture bias. Studies with simulated LHCb data documented in Ref. p3] show 
that the curvature bias is proportional to the x displacement of the detector 
planes behind the magnet. The reported constant of proportionality allows 
us to predict the displacement from the observed slopes in Fig. [5j We com- 
pute a value of approximately 1.8 mm for the x movement of the tracking 
layers behind the magnet, consistent with the actually observed difference in 
alignment constants between the scenarios with and without mass constraint. 

As explained in the introduction a curvature bias may appear in the 
scenario without the mass constraint because it corresponds to a weak mode, 
a common displacement of tracking layers that is not constrained by single 
particle trajectories. Weak modes are constrained in the LHCb alignment 
by survey information. The precision of the survey is poor compared to the 
resolution of individual tracks. However, since tracks carry no information 
on the weak modes, in the absence of mass constraints the precision of the 
survey limits the precision to which the curvature bias can be constrained. 
The 2 mm difference between the two alignments seems large compared to the 
typical uncertainty of the survey information for the detector layers behind 
the magnet, which is about 0.5 mm. However, for the results obtained here 
the vertex detector was used as a reference. The observed displacement is 
equivalent to a 0.2 mrad rotation of the vertex detector around its y axis, 
consistent with the precision of the vertex detector survey. 



5. Conclusion 

Following a recipe outlined in Ref. [B] we have developed a method to 
extend a track-based minimum \ 2 algorithm for detector alignment with 
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information from reconstructed vertices. We have presented the algebra of 
the vertex fit, including a mass constraint. We have demonstrated how such 
a vertex fit can be exploited to remove the effect of a curvature bias in the 
alignment procedure. For LHCb, invariant mass resolution is important to 
isolate rare B decays from the background and to separate decays that are 
kinematically close, such as decays of B mesons to two light hadrons. Mass 
constraints are now part of the standard LHCb alignment procedure and 
allow the LHCb collaboration to perform world-best measurements in such 
decays [T7] . 
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Appendix A. Track and vertex model for a forward detector 



Here we discuss a parametrization for the vertex fit in a forward detector 
like LHCb. We define a cartesian coordinate frame as in Section HI Track 
trajectories are locally parametrized by a 5-D 'state' vector (x,y,t x ,t y , k) at 
a given z coordinate, such that (x, y, z) is a point on the track, t x and t y are 
the local tangents dx/dz and dy/dz, respectively, and k = Q/(pc) with Q 
the charge and p the momentum. 

In the vertex fit the vertex coordinates £ are parametrized by a vector 
£ = (x v ,y v ,z v ). The momentum vector of each of the outgoing tracks i is 
parametrized by a 3-D vector $ = (t X)V> i,ty )Vt i,K V) i), where the meaning of 
parameters is the same as in the measured track state. The subscript v 
indicates that these are now parameters of the fit, not reconstructed track 
parameters. The motivation for choosing this particular parametrization is 
that it makes the measurement model nearly linear. 

We assume that the reference z coordinate of each of the track states is 
sufficiently close to the z coordinate of the vertex that the magnetic field can 
be ignored. The measurement model for track i then becomes 
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where Z{ is the position at which the measured state vector of the track is 
defined. For the corresponding derivative matrices we obtain 
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(A.3) 



Note that in the absence of a magnetic field, the magnitude of the momen- 
tum is not required to obtain the vertex position. Therefore, the fit can be 
performed more efficiently by omitting the curvature parameter altogether 
and working with a 4-D track model. The updated curvature can then be 
computed afterward by propagating the change in the remaining four track 
parameters. This approach requires more algebra, in particular when imple- 
menting the mass constraint. Consequently, we have not used it. 

For completeness we also present the expressions for the derivative matri- 
ces needed for the mass constraint. We drop the subscript v for readability 
and choose units such that c = 1. In terms of the fit parameters the four- 
vector (ki, ko t i) of outgoing track i becomes 



/ Pi t x ,i / Hi \ 
Pi ty,i I Tli 

Pi /rii 
V V m i + P 2 i J 



(A.4) 



where is the track candidate mass, Pi = Qi/K,i and ni = y/T+ t Xii 2 + t 
The derivative matrix in Eq. [21] is then given by 



y-i 



/ (1 + tl i) pi I nf -t X)i t yA pi I n\ 



-t x ,i ty :i Pijn\ (1 + t 2 x i ) pi I n\ 



\ 



~tx,i Pi I Tlj 





~ty,i Pi I 





-Qit x ,ip 2 I rii \ 

Qi ty,i Pi I Tli 

-QiPi /nj 

-QiPi I V m t + P 2 i J 



(A.5) 



Finally, taking the vertex invariant mass as the norm of the total four- vector, 
k tot \, one obtains for the derivative matrix Eq. 



m = 
constraint 



^0,tot 
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of the mass 



Hm = 



(-k Xt tot/m, -ky jtot /m, -k Xjtot /m, fc ,W m ) • ( A - 6 ) 
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